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ABSTRACT 

We explore the build-up of stellar mass in galaxies over a wide redshift range 0.4 < z < 5.0 by study- 
ing the evolution of the specific star formation rate (SSFR), defined as the star formation rate per unit 
stellar mass, as a function of stellar mass and age. Our work is based on a combined sample of ~ 9000 
galaxies from the FORS Deep Field and the GOODS-S field, providing high statistical accuracy and 
relative insensitivity against cosmic variance. As at lower redshifts, we find that lower-mass galaxies 
show higher SSFRs than higher mass galaxies, although highly obscured galaxies remain undetected 
in our sample. Furthermore, the highest mass galaxies contain the oldest stellar populations at all 
redshifts, in principle agreement with the existence of evolved, massive galaxies at 1 < z < 3. It is 
remarkable, however, that this trend continues to very high redshifts of z ~ 4. We also show that 
with increasing redshift the SSFR for massive galaxies increases by a factor of ~ 10, reaching the era 
of their formation at z ~ 2 and beyond. These findings can be interpreted as evidence for an early 
epoch of star formation in the most massive galaxies, and ongoing star- formation activity in lower 
mass galaxies. 

Subject headings: surveys — galaxies: evolution — galaxies: formation — galaxies: fundamental 
parameters — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

In recent years, there has been considerable interest in 
the relation of the stellar mass in galaxies and their star- 
formation rate (SFR), since this allows to quantify the 
contribution of the recent star formation to the build up 
of stel lar mass for different galaxy masses. iCowie et all 
(1996) used iGband luminosities and [Oil], Ha or ul- 
traviolet (UV) fluxes to investigate this connection for a 
/^-selected sample of ~ 400 galaxies at z < 1.5 and noted 
an emerging population of massive, heavily star forming 
galaxies at higher redshifts, a phenomenon they termed 
'downsizing'. Later on, the 'specific SFR' (SSFR), de- 
fined as the SFR per unit stellar mass, was used to study 
t his relation. 

iGuzman et all l)1997j) derived the SSFR for 51 com- 
pact galaxies at z < 1.4 in the HDF flanking fields 
( Willi am's et al.ll996f l finding no evi dence for an increase 
of the peak SSFR with redshift. iBrinchmann fc Ellis] 
l|2000|) studied 321 /-selected field galaxies at z < 1 
and detected a clear upper limit on the SSFR moving 
to higher SFRs with increasing redshift. They conclude 
that the most massive galaxies must have formed the 
bulk o f the ir stars before z = 1. | Perez-Gonzalez et all 
(2003) and IBrinchmann et al.l l|2004h presented detailed 
investigations of the SSFR in the local universe, while 
iFontana et al.l l)2003j) used a deep if-selected sample of 
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~ 300 galaxies in the HDF-S l|Casertano et al.M2000|) to 
trace the SSFR to z > 2 and found more evidence for 
hi gher SSFRs in the past, a result confirmed at z < 1.5 
bv lBauer _e^frlJ_j|2QQ5j)_us ing sp ectroscopic data for ~ 350 
galaxies. Uuneau et alJ l)2005[) found clear evidence for 
downsizing in a sam ple of ~ 200 galaxies at 0.8 < z < 2 
from the OD DS. andlFeulner et all ll2005h used the M U- 
NICS data llProrv et a,l.l 120011 iFeulner et alJ 12003^ to 
study the SSFR of ~ 6000 galaxies with photometric 
redshifts to z — 1.2, placing strong emphasis on the age 
of the stellar populations. They confirmed previous re- 
sults on the rise of the SSFR with redshift, but found 
in addition that the highest mass galaxies are dominated 
b y the oldest stel l ar po pulations at all redshifts. 

Ha mmer et alJ |200 5'1 obtained 15/mi fluxes for ~ 200 
z > 0.4 galaxies and estimated that 15% of all Mb < —20 
galaxies are luminous IR galaxies with SSFRs well above 
the range usually found usin g other star-formation es- 
timators. iBell et all l|2005bD investigated ~ 1700 B- 
selected galaxies at z ~ 0.7 with photometric redshifts, 
~ 25% of which could be detected at 24/im. They found 
that these galaxies typically have masses in the range 
9.5 < M„ < 11.0 and SFRs of up to < 100 M yr" 1 . 
iPerez-Gonzalez et al.l (|2005|) used about 8000 sources se- 
lected at 24/im to study the SSFR to z ~ 3, finding clear 
support for the downsizing picture. 

In the following we present results on the SSFR of 
- 9000 galaxies at 0.4 < z < 5.0 in the FORS Deep Field 
(FDF) and the GOODS-S field, reaching higher redshifts 
than previous investigations. This letter is organized as 
follows: We introduce the galaxy sample and our method 
to derive the SSFR in Sectional In Section|2|we present 
our results on the evolution of the SSFR with redshift. 
Section 0] gives a brief account on the influence of dust 
attenuation. In Sectional we discuss different evolution- 
ary paths in the SSFR-stellar mass diagram, before we 
summarize our findings in Section Throughout we as- 



2 



Georg Feulner et al. 



O 



■ 0.01 Gyr 


■ 1.50 Gyr 


■ 0.05 Gyr 


■ 2.00 Gyr 


■ 0.10 Gyr 


■ 2.50 Gyr 


■ 0.20 Gyr 




■ 0.40 Gyr 


■ 4.00 Gyr 


■ 0.60 Gyr 


■ 5.00 Gyr 


■ 0.80 Gyr 


■ 6.00 Gyr 


■ 1.00 Gyr 


■ 8.00 Gyr 




9 10 11 12 



Log M star [M sun] 

Fig. 1. — Six panels to the left: The SSFR as a function of stellar mass and redshift for the FDF and GOODS-S. Both the SFR and 
the mass are corrected for dust extinction. The solid and dashed lines correspond to SFRs of 1 Mg yr — 1 and 5 Mq yr — x , respectively. 
Objects are colored according to the age of the main component of the stellar population synthesis model fit to the photometry, ranging 
from 0.01 Gyr (purple) to 8 Gyr (red). The dot-dashed line is the SSFR required to double a galaxy's mass between each redshift epoch 
and today (assuming constant SFR); the corresponding look-back time is indicated in each panel. The error cross in each p anel gives an 
idea of the typical errors , while the dotted line roughly represents the high-mass cut-off of the local stellar mass function IDrorv et alJ 
2004; Fontana et al. 2004; Drorv ct al. 2005). Lower right-hand panel: Examples for evolutionary paths yielding a doubling of a galaxy's 
mass, through quiescent star formation, through a burst of star formation superimposed on quiescent star formation, and through a dry 
equal-mass merger. Open symbols denote the starting point, filled symbols the final state; the doubling line is drawn for a lookback time 
of 10 Gyr (z ~ 2). The arrows indicated the influence of gas consumption or loss. 



sume fi m = 0.3, Q\ — 0.7 and H = 70kms 1 Mpc . 
All magnitudes are in the Vega system. 

2. THE GALAXY SAMPLES 

The FDF ijHeidt et alJl2003() offers photometry in the 
U, B, g, R, I, 834 nm, z, J and K bands a nd is com- 
plime nted by deep spectroscopic observations ijNoll et alJ 
I2004J) . In this letter we use the /-selected sub-sample cov- 
ering the deep central part of the fiel d (~ 40 arcmin 2 ) 
as described in lGabasch et alJ l|2004a|) . containing 5557 
galaxies down to I — 26.4. 

Our X-band selected catalog for the GOODS-S field 
(Salvato et al., in prep.) is based on the publicly avail- 
able 8 2.5 x 2.5 arcmin 2 J, H, Ks VLT/ISAAC im- 
ages and contains 3297 galaxies down to K = 23.5 
in ~ 50 arcmin 2 . The U and / images are described 
in lArnouts et all (1200 1T> . the B, V, and R images in 
lsThmn^r^t^uT " l|2003|) . The data were analyzed in a 
ve ry similar way to the FDF data and already used 
in iGabasch eTaLl (|2004bf) and IDrorv et alJ p005|) . Al- 
though both samples in itself are not large in area, hav- 
ing two different lines of sight helps to overcome some 
of the effects of cosmic variance. Furthermore, we chose 
our redshift intervals large enough to further minimize 
the effect. 

Photom etric redshifts a re der ived using the method de- 
scribed in iBender et al.1 ((20011) . We estimate the SFRs 
of our galaxies from the spectral energy distribution by 
deriving the luminosity at A = 1 500 ± lOOA and conv ert- 
ing it to an SFR as described in lMadau et al.1 (|1998f) as- 



suming a Salpeter initial mass function l)Salpeterl fl955). 
Although this is an extrapolation for the lower redshift 
bins, the results agree very w ell with our work at lower 
redshifts (Feulner ct al. 2005). Stellar masses are com- 
puted from t he multi-color phot ometry using the same 
method as in IDrorv et al.1 1)2005(1 . It is described in de- 
tail and tested against spectroscopic and dy namical mass 
estimates in Drory, Bender & Hopp ((200411 . In brief, we 
derive stellar masses by fitting a grid of ste l lar po pula- 
tion synthesis models bv lBruzual fc CharlotJ l(2003() with 
a range of star formation histories (SFHs), ages, metal- 
licities and dust attenuations to the broad-band photom- 
etry. We describe the SFHs by a two-component model 
consisting of a main component with a smooth SFH 
oc exp(— t/r) and a burst contributing up to 15% in mass. 
We allow SFH timescales r £ [0.1, oo] Gyr, metallicitics 
[Fe/H] G [—0.6,0.3], ages between 0.01 Gyr and the age 
of the universe at the objects' redshift, and independent 
extinction values of Ay £ [0.0, 1.5] for the main compo- 
nent and the burst, respectively. We adopt a Salpeter ini- 
tial mass function for both components, with lower and 
upper mass cutoffs of 0.1 and 100 Mq. The SFR is cor- 
rected with the dust attenuation obtain ed for the burs t 
component using the extinction curve of lCalzettil (fl997^ . 
Note that at higher redshift the uncertainty in the mass 
estimate increases since the observed K band t hen probes 
the rest- frame blue or UV l|Drorv et alJ l2005). We have 
verified the uncertainty in the mass estimate by compar- 
ing masses at lower redshifts from simulations with and 
without the NIR bands. 
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3. THE SPECIFIC STAR FORMATION RATE 

One way to explore the contribution of star formation 
to the growth of stellar mass in galaxies of different mass 
is to st udy the redshift evolution of the specific SFR 
(SSFR; IGuzman et aTllT997t IBrinchmann fc Ellisl 1200(1) 
which is defined as the SFR per unit stellar mass. In 
Figure ^ we present the SSFR as a function of stellar 
mass and age for six different redshift bins covering the 
range 0.4 < z < 5.0. We have convinced ourselves that 
the distributions of galaxies from the FDF and GOODS- 
S are in very good agreement. 

Several effects can be observed in Fig. ^ The 
upper cut-off of the SSFR running essentially paral- 
lel to lines of constant SFR and shifting to higher 
SFRs wit h increasing redshift was already noted in ear- 
lier work llBrinchmann fc Ellisil2000T iFeulner et alJl2005t 
iBauer et al.ll20 05). This trend seems to continue to the 
highest redshifts probed by our sample: While at z ~ 0.6 
we find SFR max ~ 5 M© yr -1 , galaxies reach as much as 
SFR max ~ 100 M Q yr" 1 at z ~ 4. Note that this up- 
per envelope is partly due to a selection effect: Heavily 
dust obscure d star bursts cannot be detected in our sam- 
ple (s ee, e.g.. lHammer et al"ll2005UPerez-Gonzalez et al.l 
2005), but our conclusions still hold for galaxies not heav- 
ily affected by dust extinction (see the discussion below) . 
Furthermore, it is evident from the distribution of ages 
in this diagram, that the most massive galaxies contain 
th e oldest stellar popul ation s, as has been show n already 
in iFeulner et ail <)2005D and iDrorv et all 1200511. This is 
in ag reement with the 'downsizing' scenario ijCowie et al.l 
1996) 

4. THE ROLE OF DUST 

The influence of dust extinction on the determination 
of the SSFR is two-fold. First, heavily dust enshrouded 
objects might escape detection because too much of the 
optical light is absorbed. Secondly, objects might be de- 
tected, but their SFR (and stellar mass) might be un- 
derestimated because of the increasing dust extinction 
in the UV. We try to correct for the second effect by 
including dust attenuation in our model fitting, and cor- 
recting both the SFR and the stellar mass accordingly, 
but the correction will likely be underestimated for ex- 
tremely dusty objects. It is more complicated, of course, 
to overcome the first effect. 

In principle both sources of uncertainty could be over- 
come with observations in the thermal infrared (IR), 
where the radiation absorbed by the dust component 
is re-emitted. Note however, that this approach suffers 
from confusion and identification problems and also in- 
volves uncertainties in the conversion of the observed IR 
flux to the total IR flux, and in the unknown contribu- 
tion of dust heat ing by old stellar p opulations, see, e.g., 
the discussion in lBell et"al] l)2005b() . 

Luminous IR galaxies (LIRGs) are a well known 
population of dusty galaxies w ith very high SFRs 
and thus also SSFRs (see, e.g., lHammer et all 120051 
iPerez-Gonzalez et al.ll2f)05jl . They can be interpreted as 
galaxies experiencing a brief episode of heavy star for- 
mation triggered by mergers or gas infall. Due to their 
limited gas supply, these galaxies would spend most of 
their time in a 'normal' state with lower SSFR. There- 
fore, although our survey misses dust enshrouded star 
forming galaxies, they can be considered as intermittent 
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Fig. 2. — Average SSFRs for galaxies with stellar masses of 
log M^/Mq G [8.5,9.5] (blue), [9.5,10.5] (green) and [10.5,11.5] 
(red) and SFRs larger than 1 Mq yr — 1 as a function of z for 
FDF (open squares), GOODS-S (open diamonds) and the com- 
bined sample (filled circles). The error bar represents the error of 
the mean. All numbers are given in Table Q] The solid line indi- 
cates the doubling line of Fig. which can be used to discriminate 
quiescent and heavily star forming galaxies. 



stages in the evolution of galaxies. In particular our re- 
sults on the existence of massive evolved galaxies even at 
high redshifts remain unaffected. 

5. EVOLUTIONARY PATHS 

In order to understand the role of various evolutionary 
paths in Fig. ^ h is helpful to visualize schematically 
different ways to double a galaxy's mass as shown in the 
lower right-hand panel of Fig.^ In the following we will 
discuss these paths in more detail. 

Quiescently star-forming galaxies: A galaxy doubling 
its stellar mass by quiescent star formation at 0.5M Q yr _1 
moves along a line of constant SFR towards the lower 
right part of the diagram. Note that galaxies below the 
doubling line in Fig. ^ do not have enough time to double 
their mass until today. 

Starbursts: In contrast to a quiescent galaxy, a star- 
burst can increase its mass in a shorter time interval, 
provided it has enough gas to consume. Bursts of star 
formation may be triggered by gas inflow or galaxy in- 
teractions, and quickly move a galaxy to high SSFRs, 
where it stays for a brief period of time before it fades 
back to 'normal' SSFRs. Given the typical dusty nature 
of starburst galaxies, they might escape detection in opti- 
cal surveys during this stage. However, since these bursts 
are typically brief, the galaxies spend most of their time 
with the quiescent galaxies. Note that the exact path de- 
pends on the details of the star formation history during 
the burst phase, however, the only relevant parameter 
for the final mass is the mass of the consumed gas. 

Dry mergers: Two galaxies undergoing a dry merger 
(i.e. without interaction-induced star formation) basi- 
cally move to the right in the diagram. We illustrate this 
with equal stellar mass mergers; the stellar mass clearly 
doubles, while the final SSFR is the average of the two 
initial SSFRs. 

Note that the 'true' endpoint of the galaxies' evolution 
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fe.g. lSmail et alJl2002|) . might shift this formation epoch 

TABLE 1 t lower redshifts, but does not affect our conclusions. 

Average SSFRs of galaxies as a function of z 



(z) log SSFR log SSFR log SSFR 

(1) (2) (3) 



0.60 


0.40 ± 0.04 


-0.20 ± 0.03 


-1.11 ± 0.10 


1.00 


0.40 ± 0.02 


-0.09 ± 0.03 


-1.04 ± 0.06 


1.50 


0.41 ± 0.02 


-0.07 ± 0.02 


-1.02 ± 0.08 


2.15 


0.59 ± 0.03 


0.16 ± 0.03 


-0.56 ± 0.07 


3.00 


0.62 ± 0.03 


0.19 ± 0.02 


-0.21 ± 0.08 


4.25 


0.75 ± 0.04 


0.22 ± 0.04 


-0.20 ± 0.08 



Note. — The SSFR is given in units of Gyr 1 . 
(1) log M*/M G [8.5,9.5]. (2) logM*/M € 
[9.5,10.5]. (3) logM*/M Q e [10.5,11.5]. 



will in all three cases likely be lower than shown, since 
all three processes diminish the limited gas supply. This 
is indicated by the down-ward pointing arrows. 

In the light of these evolutionary possibilities it is ob- 
vious that the only two ways to form massive galaxies 
with old stellar populations is by highly efficient early 
star formation in massive haloes, or by dry merging of 
less massive galaxies harboring old stars. Both scenarios 
can, in principle, be distinguished by analyzing the red- 
shift dependence of the SSFR in the most massive galax- 
ies. The result is presented in Fig. [21 where we show 
the average SSFR as a function of redshift for galaxies in 
three mass intervals. At redshifts z < 2, the most mas- 
sive galaxies with log M*/Mq <E [10.5, 11.5] are in a qui- 
escent state with SSFRs not contributing significantly to 
their growth in stellar mass. However, at redshifts z > 2, 
the picture changes dramatically: The SSFR for massive 
galaxies increases by a factor of ~ 10 until we witness the 
epoch of their formation at z ~ 2 and beyond. The fact 
that we miss galaxies at z > 1 which have both high star- 
formation rate and mass, e.g. sub-mm galaxies at z ~ 2 



6. DISCUSSION AND CONCLUSIONS 

It is remarkable that the most massive galaxies show 
evidence of harboring the oldest stellar populations at all 
redshifts. Our sample shows this effect robustly out to 
very high redshifts of z ~ 4. Note that this is in agree- 
ment with the findings of massive, evolved galaxies in the 
p opulation of Extremely Red Objects (EROs) at 1 < z < 
2 llSaracco et al.ll2003[lOimatti et alJl2004HSaracco et alJ 
120051 lLonghetti et al.l[2005D an d among: the Distan t Red 
Galaxies (DRGs) at 2 < z < 3 (jLabbe et al.ll2005|) . Ap- 
parently this trend continues to even higher redshifts, in- 
dicating a very early formation epoch for the most mas- 
sive gal axies in the unive rse, favoring the 'downsizing' 
picture ijCowie et al.lll996|) . 

This important finding is evident in Fig. where we 
show the average SSFR of galaxies with different masses 
as a function of redshift. While at redshifts z < 2 the 
most massive galaxies are in a quiescent state, at red- 
shifts z > 2 the SSFR for massive galaxies increases by 
a factor of ~ 10 reaching the epoch of their formation 
at z ~ 2 and beyond. While there i s evidence for dry 
merging in the fi eld galaxy population l|Faber e t al. 2005; 
iBell et al.l2005aj) . this strong increase in the SSFR of the 
most massive galaxies suggests that at least part of this 
population was formed in a period of efficient star for- 
mation in massive haloes. 
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